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Abstract 

The next-to-leading-order (NLO) 0(1/(5^) power correction for 777 and tj'^ form factors are 

evaluated and employed to explore the r] — ij' mixing. The parameters of the two mixing angle 

I scheme are extracted from the data for form factors, two photon decay widths and radiative J/ip 

' decays. The analysis gives the result: /^^ = (1.16 ± O.OG)/,^, /»78 = (1.33 ± 0.23)/7r,^i = 
Ph. 

. —9° ± 3°, 0s = —21.3° ± 2.3°, where /rj^gj and ^^(g) are the decay constants and the mixing angles 

I for the singlet (octet) state. In addition, we arrive at a stringent range for : —10 MeV< f^, < —4 
MeV. 
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I. INTRODUCTION 



Recently, the next-to-leading order (NLO) power correction has been shown an important 
role in understanding the exclusive processes 7*7r ^ 7 and 7*7r ^ tt [Q. The method for 
calculating the NLO power correction is called the collinear expansion ||, ^, ||. This 

power expansion method is compatible with PQCD factorization which demonstrates 

the amplitude as convolutions of perturbatively calculable hard scattering amplitudes (the 
hard function) and nonperturbative hadronic wave functions (the soft function). Further- 
more, it is a Feynman diagram approach such that the partonic interpretation for the NLO 
power correction can be preserved. 

The asymptotic limit of the 717 transition form factor to be with = 93 MeV, 

is about 15% higher than the upper end of the CLEO data 0. This deviation can be 
successfully explained by NLO power correction. In addition, NLO power correction can 
well describe the low energy portion of the CLEO data. Our purpose in this paper is to 
generalize the approach for the 717 form factor to the 777 and Vj^ form factors. The first 
problem we shall face is that there are many independent degrees of freedom associated 
with the ri and rj lowest valence Fock states. The ri and Vj mesons are admixtures of 
the SU{?i)p octet and singlet states, this counts eight quantities: four wave functions and 
four related decay constants. The // and Vj mesons can also receive contributions from the 
\J{\)a anomaly to have intrinsic heavy quark and gluon contents. To reduce the number 
of independent degrees of freedom, we invoke phenomenological constraints and physical 
assumptions over the wave functions and decay constants. 

For the phenomenological constraints, we shall employ the large momentum transfer data 
for ?77 and Vj^ transition form factor |, 0, , the two-photon decay widths of the 77 
and rj mesons and the ratio -Rj/^ of the J/ip ^ 7]'-) and J/ip ^ rj^ decay widths. As for 
the physical assumptions we shall invoke the SUi^^p octet-singlet mixing scheme. In this 
mixing scheme, both 77 and r( are linear combinations of r/s and r/i, the octet and singlet 
states in the SlJi^^p representation. The mixing is controlled by the mixing angle. In the 
one mixing angle scheme, in which the octet and singlet meson states have different decay 
constants /.^g and /^^ and share a common mixing angle = Q\ = 0, the mixing angle Q 



is in range of —20° < 6 < —10° |T2|. In recent years, many evidences [|1^, |T8|] have 
indicated that a two mixing angle scheme, in which the octet and singlet mixing angles 
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and 6i can take different values, is more general than the one mixing angle scheme. The 777 



and rf'^y form factors have been investigated by either one mixing angle scheme |]T5|, [T6|, [T7 
or two mixing angle scheme [0]. We would like to make a more refined analysis for these 
form factors by including NLO power correction. 

In this paper, we shall investigate the 777 and rj''^ form factors and the mixing pattern 
oi 7] — 7]' system by employing PQCD formula with NLO power correction in the standard 
hard scattering approach. We shall employ the coUinear expansion to derive NLO (0((5~^)) 
power correction for 777 and r]''-) form factors in Sec. II. In Sec. Ill we shall analyze the high 
momentum transfer data for form factors, the two photon decay widths and the ratio -Rj/^. 
The values of the mixing parameters are determined from a analysis for the data. Sec. IV 
is devoted to conclusions. 



II. r/7 AND 7?'7 FORM FACTORS AND COLLINEAR EXPANSION 

Our strategy in calculating the power corrections to the 77 and rj' meson-photon transition 
form factor is to invoke the coUinear expansion |^, ^ ^ . For simplicity, we shall first ignore 
the meson mass effects. That is we choose the momentum of the initial state meson. Pi, 
and that of the final state photon, P25 as 

/\/f2 Juf'i 

P2' = (0,|,0^)^^7^^ (1) 

such that the virtual photon has momentum q = P2 — Pi with virtuality = —Q^ to make 
PQCD applicable. Vectors p and n are in the + and — directions in the light-cone reference 
frame and have properties = = and p ■ n = 1. Mp denotes the mass of the initial 
state meson. For the Feynman diagrams displayed in Figs. 1(a) and (b), the amplitudes are 
written as 

A{Pi, P,) = I ^Jr[H{k, P2, Q'Mk, Pi, Q')] (2) 

where the trace is taken over the color and spin indices and the meson DA Pi, Q^) has 
expression for P(= 77,77') meson 

$(fc,Pi,Q2) = / c/Se^'-^(0|g(0)g(z)|P) . (3) 
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We assign the loop momentum k for the valence antiquark and let it flow into the hard 
function. The hard function H{k, P2, Q^) contains two parton photon interaction vertices 
and one virtual internal parton propagator. The amplitude A contains leading, next-to- 
leading and higher twist contributions. The quantity twist is understood as an effective 
twist for nonlocal operators and is not exactly the same with the usual twist defined for 
local operators. The twist has different meanings for the hard and soft functions. For 
the hard function, the twist is defined as the power of the inverse of the photon virtuality 
Q, and, for the soft function, the twist represents the power of the small scale A with 
magnitude of order Aqcd- By employing collinear expansion, we can systematically separate 
the leading twist (LT) contributions from the next-to-leading twist (NLT) contributions. The 
LT contributions are from collinear loop momentum k = xp. It is therefore convenient to 
parameterize the loop momentum k into 



k^ = xp^ + ^^ n^ + k^, , (4) 
2x 



where k contains on-shell part 



and off-shell part 



^2 

kl = xp^ + -^n^ + k1 (5) 



k^s = . (6) 



In the first step, we expand the hard function H{k) with respect to k as 

dH{k) 



H{k) = H{k = k) + 



{k~kY + --- (7) 

k=k 



dk"^ 

With the help of k^ and ^5, we can factorize the loop parton propagator F{k) = 
k — ie) into its long distance part, Fiik), and short distance part (the special propagator 
defined in 0), Fs{k), which take expressions 

FLik) = ^ , Fsik) = -f^ . (8) 
k^ — It 2k ■ n — te 

The propagators Fi{k) and Fs{k) have different physical meanings. To see this, it is amount 
to consider their propagations on the light cone. The integrals of Fi{k) and Fs{k) over 
k~^ = k ■ n give 



where t] and A mean the hght-cone distances in the — direction. It is obvious that Fs{k) 
is not propagating on the hght cone. This means that Fs{k) should be included into the 
hard function. By dimensional counting, Fs{k) is of order 0{Q^^). Therefore, including one 
Fs{k) into the hard function then increases one twist order for the hard function. 

There are different effects as -Fl(^) and Fs{k) act on the spin structures of hard function, 
the terms proportional to ^ or fi. As Fi{k) acts on its collinear part vanishes and non- 
collinear parts are retained 

FL{k)i> = -F{k){k - kr{t^^)Fsm (10) 

where minus sign comes from the anti-particle propagator. The vertex i'-ja and short distance 
propagator Fs{k) are then absorbed into the hard function. The factor {k — k)°' is included 
into the soft function to become a coordinate derivative on the quark fields. As -Fl(^) acts 
on fi, its collinear part contributes to leading order. The short distance propagator Fs{k) 
only serves to introduce the interaction term qjiq for ^ vertex, where A°' denote the gluon 
fields. The total effects of Fiik) and Fsik) acting on |^ are to include one i-^a and one Fsik) 
into the hard function and to absorb the factor [k — A;)" and gauge fields A"' into the soft 
function to become a covariant derivative, -D" = id"' — gA" with g the strong coupling. 

The contributions from the second term of Eq. (|^ and from Figs. l.(c) and (d) are of 
twist-6 or higher twist and will not be considered in below discussions. The reason is that 
the possible non-vanishing components of in dH{k)/dk°' are a = -|- or — , but both vanish 
as dH{k) /dk" contract with {k — k)" or {0\qA"'q\P). We substitute the first term of Eq. 
into the integral with the soft function and apply the identity 

J dx5{x — k ■ n) = 1 (11) 

to convert the loop momentum integral into the fraction variable integral. The amplitude 
then becomes, approximately, 

A{Pi,P2) ^ j dxTT[H{x)<^{x)] (12) 

which contains LT and NLT contributions. The meson DA $(a;) has the expression 

/°° fix 
^e^^\0\q{0)q{Xn)\P) . (13) 

We now discuss how to separate the LT from the NLT contributions for amplitude A. Due to 
the fact that the final state photon is real and has transverse polarization, the hard function 
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H{x) can have spin structures: 7±7i7_L, 7±|^|4 7_l|4|^ and 7_l|^7_l where ■j± = 7° with a = 1,2. 
The first spin structure leads to LT contribution, while the second and third ones result 
in the NLT contributions. The last spin structure would lead to next-next-to-leading twist 
contribution and will not be considered below. To calculate the NLT contributions, we need 
to apply Eq. (|I^) to extract the contributions from non-coUinear loop momentum. As a 
result, we get the amplitude up to NLT as 

J dxTr[H{x)^{x)] + y"dxTr[i7«(x)<,$"'(x)] (14) 

where the first term of the right hand side of Eq. (|l^) comes from the Feynman diagrams 
shown in Fig. 1(a) and (b) and the second term from those diagrams shown in Fig. 2. The 
tensor w^, is defined as w^' = g^' ~ P^'na'- The NLT hard function Ha{x) is defined as 

H.ix) = {t^^)^H{x) + /f(x)(z7.)^^^^ (15) 

and the NLT meson DA $°(x) has expression 

$"(x) = dxi — ^e'("i-")''e'"^(0|g(0)D°(r7n)g(An)|P) . (16) 
io Jo 27r Jo 27r 

The factorization of momentum integral is finished. To complete the factorization, we still 
need to perform the factorizations of the color and spin indices. To separate the color 
indices, we take the convention that the color factors of the hard function are extracted and 
absorbed into the soft function. As for the spin indices, we employ the expansion of the soft 
functions into their spin components 

r 

$" = 5^r0r^", (17) 

r 

where F denotes gamma matrix and 0^^'"^ is the related spin component of the distribution 
amplitude. For a given order of l/Q^, we choose the component 0^*^'"^ with lowest twist. 
The determination of the lowest twist 0'"'^'") can be done as follows. Firstly, we notice that 
the tensor structure of 0^^'"^ can be expressed in terms of p, n, d°f = g"^ — p"'n^ and 
= e'^^^'^p^nx- The vectors p and n have dimensions [p] = 1 and [n] = — 1 with respect to 
the hard scale Q. Secondly, note that the matrix element for the soft function is written as 

$ ~ {0\qq\P) , 

~ (0|gD"g|P) . (18) 



By the above facts, we can derive a power counting rule as follows. Consider the ■ "s 

has the fermion index F and the boson index B. The fermion index F arise from the spin 
index factorization for 2F fermion lines connecting the soft function and the hard function 
and the boson index B denotes the nj^ power of momenta in previous collinear expansion 
and the gluon lines as B = tid + no- We may write 



i 

where A denotes a small scale associated with DA. Spin polarizers Cj denote the combination 
of vectors p'^, and 7^. Variable represents the twist of DA 0*. The restrictions over 
polarizers g^'i^ ^'f 

_ Q (20) 

which are due to the fact that polarizers Cj are always projected by w^'- The dimension of 
^fii---tip;ai---aB jg (determined by dimensional analysis 

d{(j)) = 3F + B - I (21) 

By equating the dimensions of both sides of Eq.(^), one can derive the minimum of 

^min^2F + 5 + i[l-(-l)^] . (22) 

It is obvious from Eq.(p2|) that there are only finite numbers of fermion lines, gluon lines 
and derivatives contributes to a given power of l/Q^. 

We now demonstrate that the collinear expansion is compatible with the conventional 
approach for proving the PQCD factorization at one loop order of radiative correction. To 
show this, we consider the radiative correction for H{x) as displayed in Fig. 3(a). If the 
radiative gluon in Fig. 3(a) is collinear with momentum /° = = {Q,X'^/Q,X) , 

where X <^ Q, the lower virtual antiquark has the momentum {k — I) {^Q, IQ-, ^) with 
k = xPi. It is obvious that the virtual antiquark in the collinear region behaves similarly 
to the loop antiquark in the tree amplitude. The collinear expansion for Fig. 3(a) in the 
collinear region is the same as the expansion for the tree diagram Fig. 1(a) in the leading 
configuration. To demonstrate this, we write the integrand for Fig. 3(a) as 



ha = I (2^)4 p '' n^pFihyFihh'Fihhx] 

^jAHsamUl)] (23) 



where we have defined h — k — l,l2 — (Pi — k + l), Is — P2 — k + I, 
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^3a{l) = -^^-^bpF{h)F{k)^,] , (24) 



and F{li) — We, firstly, expand Hsa{l) 

Hsail) = Hsail = + - 0' (25) 

with i = {x — ^)Pi. Repeating the same considerations for the expansion of the tree amph- 
tude, we can recast I^a into 



'3a 



where 



= IW)I (|]4'="«"'"''"[^^^7^''('0'°'f ((2)7j . (27) 

It is obvious that both $3a(C) ^-nd ^s'aiO coUinear divergent for coUinear I. We introduce 
corresponding soft functions $q^^ and ^^^^ to absorbe $3a(0 ^3o(0- The corresponding 
tree level hard functions are denoted as H^^^ and h[^\ respectively. If the radiative gluons 
in Fig. 3(a) are soft, i.e. the gluons have momentum / = (/+, , l±) ~ (A, A, A), there are no 
effects on power expansion. This is because the eikonal approximation up to 0{1/Q^) can 
be applied to factorize the soft gluons from the valence quark propagator. The other two 
particle reducible diagrams Figs. 3(b) and (c) can be dealt with similarly. It is noted that 
the double logarithms in Figs. 3(a) to (c) arising from the mixing contributions from the soft 
and coUinear divergences are cancelled each other. In light-cone gauge n ■ A = 0, Figs. 3(d) 
and (c) in coUinear region are more suppressed than Figs. 3(a) to (c) in coUinear region by, 
at least, After subtracting the soft contributions (the soft and coUinear divergences) 

from the one loop radiative correction diagrams Figs. 3(a) to (c), we can obtain the one 
loop corrected hard functions (LO and NLO) H^^^ and h[^\ separately. The analysis for 
the radiative corrections to Ha{x) is simple, since it only involves radiative correction and 
has no need to consider coUinear expansion. The diagrams for the radiative corrections to 
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Ha{x) are shown in Fig. 4. As a result, up to first order in radiative and power corrections, 
we can arrive at the factorized amphtudes as 

+ ^ (4'^ + Hi,'^) ® (<|.(°) + + (iff) + Hi'^) ($S°) + $W) (28) 

where the superscript indices i{i = 0, 1) denote the order of radiative correction and the 
subscript indices j{j = 0, 1) mean the order of power correction. The notation (8> represents 
the convolution integral and the trace over the color and spin indices. To prove PQCD 
factorization, we need to generalize the one loop factorization to arbitrary orders. It can be 
done straightforwardly p2 |. 



For convenience, we may write the amplitude as 

A(7*P 7) = -teh^c,fsxP^Pie^Fp^{Q^) , (29) 

where denotes the polarization vector of the final state photon. The form factors are 
expressed in terms of the octet and singlet components 

Fp,{Q') = 4'FpaiQ') (30) 

j=8,l 

where the expansion coefficients ap\P = rj,rj' ,i = 8,1 depend on the mixing scheme (see 
next section). Due to the t] — t]' mixing, we take the octet and singlet states as the basis for 
our investigation of the form factors. The superscript i = 8 or 1 denote the contributions 
from octet or singlet current (see below notation). The leading order of Fp.^{Q'^) is calculated 
from Figs. 1(a) and (b) and takes expression 

(Q^) = 4a ^ , (31) 

^ Jo Q^x{l-x) 

where the charge factors are defined as Cs = (e^ + — 2e^) / \/6 and Ci = (e^ + + e^)/ -\/3. 

The NLO of Fp^^{Q^) is evaluated from Fig.2 

^yw-)^-iaa/.. l^-''-'^^e^;-^-" . (32) 

We have taken the symmetry between the exchange of x ^ (1 — x) for 0p.(x), Gp.{x) and 
Gp.{x). The relevant DAs are expressed explicitly as follows 

1 poo J\ 

<PpXx) = -^- ^e^^^(0|J^(0,An)|P,(Pi)) (33) 

GpA^) = fdx, r^^e^''(^-^)e^^^(0|j;,(0,r^n,An)|P,(PO) , (34) 

8 Jo Jo (^tt) (27r) 

GUx) = -^rff ^'rfxi^°'^^e^^(--^)e^^^(0|j;,5(0,r^n,An)|PKPi)) . (35) 



where nonlocal currents are defined as 
1 



j|(0,An) 
4(0, An) 

Jl/^bi'^^Vn, An) 



1 

1 

7! 
1 

1 

1 

7! 



[n(0)75|4M(An) + (i(0)75|4(i(An) - 2s(0)75|4s(An)] , 
['u(0)75|4ti(An) + J(0)757i(i(An) + s(0)75|ls(An)] , 
['u(0)7Q,D/3(?7n)-u(An) + (J(0)7Q,i5/3(r7n)(i(An) — 2s{0)jaDf3{rin)s{\n)] , 
[u{0)'yaDi3{rin)u{Xn) + d{Q)'^aDf3{'r]n)d{\n) + s{Q)'^aDi3{r]n)s{\n)] , 
[«(0)757a-D/3(r?n)M(An) + d{0)-f5-faDf3{r]n)d{\n) - 2s{0)'j5'jaDp{r]n)s{Xn)] , 
[u(0)757„L'^(?7n)u(An) + d{0)'j5'jaDp{r]n)d{\n) + s(0)757„L'/3(r7n)s(An)(]36) 



Due to the factor l — 2x for Gp., Gp^ become dominate. The normalizations of (/)p and Gp are 
determined from the leptonic weak decay and the axial anomaly for Pj meson, respectively. 
This is similar to the pion case M. 



III. THE MIXING SCHEMES 

We employ the SU (3) p octet and singlet states to describe the r] — rj' system. The rj 
and rj' meson states can be described by means of the octet and singlet states |?78) and |?7i) 
through the one mixing angle scheme 

|?7) = cos6'|?78) — sin6'|?7i) 

\ri') = sin6'|?78) + cos6'|?7i) (37) 

where the mixing angle 6 controls the relative strength. With the mixing, the rj'j and r/'j 
form factors take expressions 

F^^(Q') = cos eFr,,y{Q^) - sine F^,^{Q^) , 

F„'7(Q') = sin 9 F^,^{Q^) + cos 9 F^,^{Q^) . (38) 

To proceed, we also assume that the octet and singlet DAs take asymptotical form. There- 
fore, we have (f)rjX^) = 3/^.x(l — x)/^/2 and G^i(x) = 3v^7r^/^.a;(l — x). The form factors 
F^.^{Q^) (i=8,l) are simplified by substituting 4>ni{x) and Gj^.{x) 



4Ci / dx 



x(l — x) 

r- fr, ^TT^fn 

6V2a0[i-^]. 



G„Ax] 



(39) 
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Using Eq. ( pQ] ) into Eq. (0), we may derive the coefficients Up in Eq. (pOl). 

To compare the form factors with the data, we extrapolate the form factors to all orders 



(40) 



This formula gives a theoretical support to the approach using the interpolating formula for 



the 777 and r]''y form factors |T8 

The decay constants /^., z = 8, 1 and the mixing angle 6 will be determined by a least 
fit to the transition form factor data above 1 GeV^ and the two photon decay widths 

T[ri 77] = (0.46 ± 0.04) keV , T[ri' 77] = (4.28 ± 0.19) keV . (41) 

The decay rates have theoretical expressions 

9a' 



T[ri^ 77] 



327r3 



f 1 fS _ f8 fl 
J n' J Tj J r)' J T] 

f 1 f8 _ f8 fl 

J Tj' J Tj J r}' J T] 



n 2 



(42) 



where decay constants are defined as 

= fm cos^ , fl = -U sine , 

= ^sin^ , fl, = f^^ cose. (43) 

The fit results are shown as Fit I in Figs. 5 and 6 and in Table. I. It is seen that the Fit 
I is in good agreement with the data for the form factors. To test the fit parameters, we 
employ the ratio of the decay rates for J/ ip into r/'j and 777 

It is usually assumed that the radiative J/ip v{v')l decays are dominated by nonpertur- 
bative gluon matrix elements {0\GG\r]') and (OlGGIr^) such that the ratio takes expression 

(45) 



with pp = Mj/^{1 — Mp/Mj^^)/2 being the three momentum of the P-meson. Our fit result 
is close to that one obtained from chiral perturbation theory (^PT), except the octet decay 
constant /^^ = A < 1-28/^, 



XPT : ^ -20° ~ -10° ; = 1.28^ ; ^ 1.1/^ . 



(46) 
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The octet decay constant /^^g is calculable by xPT up to one loop approximation 

M'^ M'^ 
(47r/^)^ (47r/^)^ (47r/^)^ 

It is noted that the predict results for r(?7 — > 77), V{rj' 77) and Rj/^p shown in Table I 
are close to the experimental values within Icr accuracy. 

Recently, it has been proposed |1^, |T8[ that \ri) and \7]') can mix through a two mixing 



angle scheme as 

\rj) = cos9s\r]s) - sm9i\r]i) , 

\r]') = sin 6*8 1 r/s) + cos 6*11771) , (48) 

where 6i denote the mixing angles. Using this mixing scheme, the 177 and ri''j form factors 
can be expressed, correspondingly, as 

F^^(g') = cos^8i^r,87(Q') -sin^iF^,^(g') , 

F,,^(Q2) = sin 9sFr„^{Q^) + cos 9^F,^,^{Q^) . (49) 

The form factors Fj^.^{Q'^) are the same as those in the one mixing angle scheme. We first 
change the values of /^^ and 9i to fit the data. The fit result is shown as Fit II in Figs. 5 
and 6 and Table. I. From Fit I and Fit II in Table I, it is found that the two mixing angle 
scheme is better in accuracy than the one mixing angle scheme by 100%. This is close to 
the investigations [|l^, 0, [l^ . 



A large value of = cos 9ifr,^ [psj which is responsible for the intrinsic charm content 
of the r]' meson, has been proposed to resolve the large branching ratios Br (5 v'K) and 
Br(i? — i> Xs7). We may explore this within our approach by adding intrinsic charm content 
into our formalism. The effects of the intrinsic charm content are similar to that of the singlet 
component. That means one can replace /r,^ with /^^ = f^,/ cos 6*1, the decay constant for 
intrinsic charm for the corresponding singlet terms. That is the part of the form factors 
from the intrinsic charm taking expression 



dx 



xil — x) 



(50) 



where Cc = 2/3 and the related DAs are 4>r)A^) = 3/^^x(l — x)/\/2 and Gr^^x) = 
3v^vr^/^^a;(l — x). The effects from the large value of the charm quark mass has been 
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absorbed into the twist-4 DA (^^^(x)!]!!. After including the contribution of the intrinsic 
charm, the rj'j and r]''-/ form factors then become 

Fr,^{Q^) = cosesFr,,^{Q^)-smei{Fr„^{Q^) + F^^^{Q^)) , 
Fr,,,{Q^) = sin^8F,3^(g2) + cos^i(F,,^(g2) + F,^^(g2)) . (51) 

As the case of the octet and singlet form factors, the extrapolation of F^^^(g^)) to all orders 
is implied. By observing Eq. (|5TD, form factor F^'^(g^) has a larger dependence of /^^ than 
Frj-yiQ"^)- We make a least fit to the form factor data to determine possible values of fri^ 
by keeping other parameters fixed. From Table II, one may see that including the intrinsic 
charm content can indeed improve the accuracy. This shows that our formalism is consistent 
in perturbation theory that the higher Fock state can be reasonably added in. As shown 
below that the allowed value for /^^ is less than In literature 0, ^ |2^, /^^ is 



proposed in the range —140 MeV < /^^ < 15 MeV. To test this, we plot in Fig. 7 the 
distribution for each set of mixing parameters list in Table II over a wide range of —140 
MeV < /^^ < 140 MeV. It is seen that the range of /^^: -10 MeV < /^^ < -4 MeV are 
allowed by the data. Because the value of cos 6*1 is close to unity, is almost equal to /^^. 

From the above analysis, one may see that combining the high energy data and the low 
energy experiment can result in constraints on the mixing parameters in a very efficient way. 
We can give a general analysis for the two mixing angle scheme. We shall investigate the 
distributions of the mixing parameters. The procedure of analysis is as followed. We firstly 
separate the data into two groups. The data for the form factors, the two photon decay 
rates and the ratio for J/tp radiative decays are denoted as set I while the latter two data 
are chosen as set II. We then determine the least values for sets I and II, respectively. 
The reason for separating the data into set I and set II is that neither set I nor set II can 
completely constraint the parameters. The parameters locate in la accuracy of the set I still 
have large uncertainties and require further restrictions, which can be obtained by the data 
set II. To be more explicit, we plot the allowable regions for the mixing parameters within 
la error with respect to those values associated with the minimal points. As shown in 
Figs. 8 and 9, both allowable regions for data set I and II are large while their intersections 
are rather restricted. The reason for this fact is easily understood. Within the data set I, the 
experimental errors are shared for high and low energy data. Because the distribution 
can only measure the correlations between the mixing parameters, the hope for constraining 
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each parameters in a independent way can not be obtained and only partial restrictions 
over the correlations of the parameters can be derived. This can be seen from Figs. 8, in 
which the /§ parameter is not contrained in a reasonable way. To compensate this flaw, we 
note that the distribution of the data set II can intersect with that distribution of 
the data set I. The intersections between two distributions can give better constraints 
on parameters than the data set I or II. One should note that although the data set II is a 
subset of the data set I, the distribution of the data set II is not necessary to be a subset of 
the distribution of the data set I, since the effects of correlations between the parameters 
are different for two data sets. Prom the overlapped regions of the data set I and II in la- 
error, we may extract from Figs. 8 and 9 the allowable regions for the mixing parameters. 
In Fig. 8, we plot the possible allowable ranges for fj^^ and /^^. It can be observed from 
Fig. 8 that the overlapped region for /^g and /^j is quite stringent: 1.1 < fna/ < 1-56 and 
1-1 < /tji/Ztt < 1-22. Figure 9 shows the allowable region for 9i and 9^ from both data set I 
and set II. The overlapped region indicates that: -12° < < -6° and -23.6° < ^8 < -19°. 

Combining Figs. 8 and 9, we may derive the — > oo limits of the scaled 777 and rj'^ 
form factors 

Q'Fr,^{Q^)\Q2_^^ = -^[/,3 COS ^8 - 2^/2^ sin^i] = (189 ± 46)MeV , 

Q^Fr,>^{Q^) |q2_^ = -^[^3 sin ^8 + 2^2^ cos ^1] = (295 ± 35)MeV . (52) 

The error in the scaled rj'-f form factor being larger than that of the scaled 7]''-/ form factor 
is due to the fact that the errors in the data for 777 and ri'^y form factors are shared in our 
analysis. This is consistent with the concept of rj — rj' mixing. 



IV. CONCLUSIONS 

We have shown that the coUinear expansion for j*rj{rj') 7 can be systematically 

performed in compatibility with PQCD factorization. The 0{Q^'^) power corrections for 
F-q-iiQ"^) s^iid F^i^i^Q"^) have been evaluated. The magnitudes of NLO power corrections are 
determined. 

We have made a general analysis for the allowed values for the mixing parameters by 
combining the high and low energy data. Except of /^g, the other three parameters /^j, ^8(1) 
can be constrained in a reasonable region. The large error for the fit /^g is mainly from 
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the experimental error and can be improved by future experiment with higher accuracy. At 
present, we would invoke the chiral perturbation theory calculation for /^g. 

We have also shown that the intrinsic charm content of rj meson has little contributions. 
Any sizable contribution from the intrinsic charm would lead to large value as shown in 
Fig. 7. Of course, there require further investigations for this point, if more accurate form 
factors data are available. 

So far, we have not considered the finite meson mass effects for the 777 and rj'y form factors. 
At first glance, the meson mass effects can not be safely ignored. Because the values of r\ 
and r( meson masses are large as = 547 MeV and M^/ = 958 MeV, power corrections 
from the mass effects of order 0{M'j>/Q'^) and 0{Mj,K^ /Q'^) are also important. The type 
of corrections 0{Mp/Q'^) are kinematics. This is similar to the Nachtmann's correction 
for deep inelastic scattering and can be negligible. The second type of corrections 
0(MpA^/(5^) are dynamics and can be argued that they are of twist-6, at least. This is 
because the associated spin projector is ^75 (cf. the leading spin projector 1^75), which will 
introduce two additional Fs propagators into the related hard function. As a result, the 
dynamical type meson mass corrections are of order 0{MpA'^/Q^). Of course, a complete 
analysis for these meson mass effects is important. 
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frfu 1 fir 


fr)i / fn 


^8 


Oi 


r(7? ^ 27)[keV] 


r(r?' 27)[keV] 




x7dof 


Fit I 


0.99 


1.08 


-16.4° 


-16.4° 


0.49 


4.47 


5.6 


64/31 


Fit II 


1.32 


1.16 


-22.3° 


-9.1° 


0.50 


4.34 


4.4 


32/31 



TABLE I: The results of the fit to the rjj and rj'j transition form factors and the two photon 
widths within the one and two mixing angle schemes. 





fijs 1 /tt 


/tji/Ztt 


/.c[MeV] 


^8 


Qi 


x'/dof 


Fit I 


0.99 


1.08 





-16.4° 


-16.4° 


61/29 


Fit I 


0.99 


1.08 


-8 


-16.4° 


-16.4° 


34/29 


Fit II 


1.32 


1.16 





-22.3° 


-9.1° 


31/29 


Fit II 


1.32 


I.IG 


-5.6 


-22.3° 


-9.1° 


19/29 



TABLE II: The results of the fit to the 777 and rj^ transition form factors within the one and 
two mixing angle scheme with or without intrinsic charm content. 
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(a) (b) (c) (d) 



FIG. 1: The leading order diagrams for 7*77(77') 7. The cross symbol means the vertex of the 
virtual photon. 




(a) (b) (c) (d) 



FIG. 2: The next-to-leading twist (NLT) diagrams for 7*77(77') — 7. The propagator with one bar 
means the special propagator. 




(a) (b) (c) (d) (e) 

FIG. 3: The one loop diagrams for the LT amplitude of process 7*77(77') — ^ 7. 




(a) (b) (c) (d) 




(e) (f) (g) (h) 



FIG. 4: The one loop diagrams for the NLT amplitude of process 7*77(77') —>■ 7. 
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G^[GeV2] 



FIG. 5: The results of the fit to the r/7 transition form factor within the one mixing angle 
scheme (the solid line) and the two mixing angle scheme (the dashed line). The data point are 
taken from 0, |, |, |lO|, 
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FIG. 6: The results of the fit to the r/'7 transition form factor within the one mixing angle 
scheme (the solid line) and the two mixing angle scheme (the dashed line). The data point are 
taken from 0, |, |, |lO|, 
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FIG. 7: The plot of the distribution vs. The values are evaluated for the data for forms 
factors by employing the sets of parameters list in Table I. 



FIG. 8: The plot of /^^ versus fr/g from the analysis for the data set I and II. The regions 
denoted as I and II represent the allowable values for /^^ and /^g within la error for corresponding 
data set I and II. 



FIG. 9: The plot of 9i versus 08 from the analysis for the data set I and II. The regions denoted 
as I and II represent the allowable values for 9i and Og within la error for corresponding data set 
I and II. 
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